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Abstract. The dependence of Tc for MgB2 on purely hydrostatic or nearly hydro- 
static pressure has been determined to 23 GPa for single-crystalline and to 32 GPa 
for polycrystalline samples, and found to be in good agreement. decreases from 
39 K at ambient pressure to 15 K at 32 GPa with an initial slope dT^/dP ~ -1.11(2) 
K/GPa. Evidence is presented that the differing values of dTc/dP reported in the 
literature result primarily from shear-stress effects in nonhydrostatic pressure media 
and not differences in the samples. Although comparison of these results with the- 
ory supports phonon-mediated superconductivity, a critical test of theory must await 
volume-dependent calculations based on the solution of the anisotropic Eliashberg 
equations. 



1 Introduction 

In spite of considerable efforts during the 20 months since the discovery of supercon- 
ductivity in MgB2 at 39 K |jl|, the search for a further binary compound with a higher 
value of Tc has yet to bear fruit In contrast, the high-Tc oxide Y-Ba-Cu-0 (Tc ~ 
92 K) was synthesized only a few weeks after the landmark discovery of supercon- 
ductivity near 35 K in La-Ba-Cu-0 0]. Progress with the binary compounds may 
have to wait until we first reach a clear understanding of MgB2's extraordinary super- 
conducting and normal-state properties, most of which are highly anisotropic. The 
compressibility, for example, is 64% larger along the c axis than along the a direction 
p. Under c-axis compression the electronic density of states N{Ef) is predicted to 
decrease much more rapidly (~ 0.6%/GPa) than under hydrostatic compression (~ 
0.1%/GPa) P|, the a band shifting downward and the hole number in the a band 
decreasing [0. The anisotropy in the upper critical field Hc2 increases strongly with 
decreasing temperature H, p[, approaching a value near 7 at K [|]. A large number 
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of experiments point to the existence of at least two superconducting energy gaps |TC 
which are predicted [jll], |l^, |13] to open up on different parts of the anisotropic Fermi 
surface. GonneUi et a/.|14| have very recently provided direct evidence for two-band 
superconductivity in MgB2 where the temperature-dependence of the distinct gaps 
originating from the a and tt bands were measured using point-contact spectroscopy 
on single crystals. 

The existence of multiple gaps and the highly anisotropic electronic and lattice- 



vibration properties of MgB2 have recently prompted Choi et al. [12] to attempt the 



solution of the fully anisotropic Eliashberg equations. A comparison of their calcu- 



lation with experiment has yielded promising results |12, 13|. A critical test of their 



approach would be provided by extending their calculation to reduced lattice param- 
eters, allowing a direct comparison with the results of high pressure experiments. 

Following the discovery of a new superconductor, high pressure studies are often 
among the first to be carried out. The reason for this is that the magnitude and sign 
of dTc/ dP help guide the materials scientist how to best modify the superconductor 
to raise T^. at ambient pressure, a case in point being the discovery of Y-Ba-Cu-0 
10. In addition, the high-pressure technique, sometimes in conjunction with high 
temperatures, is also invaluable to: (1) create new superconductors (e.g. sulfur 
||T5| or oxygen [|16] which metallize at Mbar pressures); (2) vary the properties of 
known superconductors (e.g. polycrystalline MgB2 [|l^ or the synthesis of MgB2 
single crystals [0, |19|, |2^); (3) induce structural phase transitions (e.g. in Ba, As, 
Bi, Sn, Ga and Tl [^); and, finally, (4) vary the lattice parameters to help identify 
the pairing mechanism and critically test theoretical models. Unfortunately, this 
final application has been traditionally underutilized by theorists. 

Hydrostatic or uniaxial pressure experiments on single crystals determine the de- 
pendence of a given property solely on the lattice parameters. Most high-pressure 
experiments, however, are carried out under nonhydrostatic conditions which subject 
the sample to a distribution of unspecified shear stresses. These shear stresses may 
be large enough to plastically deform the sample, resulting in lasting distortions and 
lattice defects. The pressure dependence of Tc may thus depend on whether the 
pressure medium is hydrostatic or not, particularly in elastically anisotropic materi- 
als like quasi ID and 2D organic superconductors or high-Tc oxides |2^. Studies 
on single crystals are of particular value since strain effects from grain boundaries in 
polycrystalline materials are avoided. 

Several studies of the dependence of Tc on pressure for polycrystalline MgB2 were 
carried out shortly after the discovery of its superconductivity ^ |2^, ^ . Using 
solid steatite pressure medium to 20 GPa, Monteverde et al. [0] reported that 
decreased under pressure at different initial rates (-0.35 to -0.8 K/GPa) for each of 
the four samples studied. On the other hand, in an experiment in fluid Fluorinert 
to ~ 1.7 GPa, Lorenz et al. [^| and Saito et al. [^| found dTc/dP ~ -1.6 K/GPa 
and -1.9 K/GPa, respectively. The first truely hydrostatic measurement of Tc{P) 
was carried out by our group to 0.7 GPa using He gas on an isotopically pure (^^B) 
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sample, revealing that decreases reversibly under hydrostatic pressure at the rate 
dTc/dP ^ -1.11(2) K/GPa ||2^; later work on isotopically pure and ^^B samples 
yielded dependences between -1.09(4) K/GPa and -1.12(3) K/GPa [gg. In further 
He-gas studies, Lorenz et al. pOf and Schlachter et al. |0 obtained -1.07 K/GPa and 
-1.13 K/GPa, respectively. Experiments utilizing fluid pressure media were carried 
out by Razavi et al. (-1.18(6) K/GPa) 0, Choi et al. (-1.36 K/GPa) [||, and 



Kazakov et al. (-1.5 K/GPa) [BH]. The flrst high-pressure measurements on a single 



crystal were carried out by Masui et al. |34| using Fluorinert, yielding dTc/dP ~ -2.0 
K/GPa. In the above experiments the values of dTd dP reported are seen to differ by 
more than a factor of two. Experiments to much higher pressures have been carried 
out by several groups and will be discussed below. The results of all known dTd dP 
measurements on MgB2 are summarized in the Table. 

In this paper we provide evidence that the variation in the reported Tc{P) depen- 
dences for MgB2 is primarily a result of shear stresses exerted by the solidifled pressure 
media on the sample. From recent measurements on high quality single-crystalline 
and polycrystalline samples we conclude that the initial dependence of the transition 
temperature on purely hydrostatic pressure is given by dTc/dP ~ —1.11(2) K/GPa 
Under nearly hydrostatic (dense He) pressures to 32 GPa [V/Vo — 0.855), we flnd 
Tc to decrease monotonically from ~ 39 K to 15 K. Although these results appear 
consistent with phonon- mediated superconductivity in MgB2, a quantitative check 
must await comparison with calculations based on the fully anisotropic Eliashberg 
equations. 



2 Experimental methods 

A wide variety of high-pressure techniques with many different pressure media have 



been used to study the dependence of Tc on pressure All techniques discussed 

below have been used at some time by our group. 

Perhaps the most widely used technique in the pressure range 1-2 GPa is a piston- 
cylinder cell in which two pistons compress a Teflon bucket containing the sample 
immersed in a fluid pressure medium such as Fluorinert, n-pentane/isopentane, or 
silicon oil. For studies to higher pressure a diamond-anvil-cell (DAG) may be used 
with a 4:1 methanol: ethanol mixture as pressure medium since this mixture remains 
fluid to approximately 10 GPa at ambient temperature. All these fluid pressure media 
freeze upon cooling and subject the immersed sample to shear stresses of varying 
strength which depend on the pressure medium used, the details of the pressure 
technique, and the rate at which the pressure cell is cooled, among other factors. 

The only pressure medium that remains fluid under high pressure (P < 0.4 GPa) 
near 40 K, where MgB2 superconducts, is helium and thus only high-pressure ex- 
periments on MgB2 in liquid helium are able to determine the dependence of Tc 
on purely hydrostatic pressure. At pressures above 0.5 GPa, helium is frozen be- 
low 39 K, but the pressure is still very nearly hydrostatic since solid helium is the 



3 



softest solid known; in addition, if the proper cooling procedure is followed, a sin- 
gle crystal of helium can be grown around the sample, resulting in negligible shear 
stresses which permit even de Haas van Alphen measurements on single crystals at 



very low temperatures ||36|. High-pressure studies in liquid helium are thus the 



measurements-of-choice for exacting quantitative studies of the properties of solids 
under high-pressure conditions. 

At the other end of the spectrum, in some high pressure techniques solid pressure 



media (e.g.steatite, NaCl or NaF) are used With a solid pressure medium. 



the application of pressure may subject the sample to relatively large shear stresses, 
possibly strong enough to plastically deform the sample or even crush it, thus intro- 
ducing a large number of lattice defects. The magnitude and direction of the shear 
stresses depend, among other things, on the pressure medium used, the temperature 
at which the pressure is changed, the pressure range, and whether the ring contain- 
ing the solid pressure medium is supported by a "belt" or not. With this so-called 
"quasihydrostatic" technique Tc may not be a reversible function of pressure. For 
these reasons high-pressure techniques using solid pressure media should be avoided 
in quantitative studies, particularly when studying elastically anisotropic materials 
such as the high-Tc oxides |2^, organic superconductors [|^, and MgB2. 

In the present experiments, dense helium is used as pressure medium both in a He- 
gas pressure cell to 1 GPa and in a DAC to over 30 GPa. The pressure in the cell is 
measured at temperatures within a few degrees of the transition temperature Tc, the 
superconducting transition is detected in a sensitive ac susceptibility measurement 
Details of these pressure techniques have been published elsewhere . Isotopically 



pure polycrystalline MgB2 samples were synthesized at the Argonne National Labs 
1^, Single crystals were grown at the Superconductivity Research Laboratory 

(ISTEC) in Tokyo, Japan; following the crystal growth at elevated temperatures, 
"type A" ("type B") crystals were quenched (slow cooled) |T^. The dimensions of 
the crystals studied in the He-gas cell were approximately 0.15 x 0.4 x 0.03 mm^ and 
in the DAC ~ 0.09 x 0.09 x 0.03 mml 



3 Results of experiment 

In Figs. 1-3 we show the results of our measurements of the pressure dependence 
of Tc on MgB2 single crystals. The midpoint of the superconducting transition in 
the real part of the ac susceptibility x' is used to define Tc [^, as seen in Fig. 1. 



The two "type B" crystals have higher values of (~ 38.2 K) than the "type A" 
crystal (~ 37.2 K). For both "type A" and "type B" crystals T^ is seen in Fig. 2 
to decrease reversibly and linearly with hydrostatic (He-gas) pressure. As with the 
polycrystalline samples, the same Tc{P) dependence is obtained whether the pressure 
is changed at low temperatures or at room temperature. As seen in the Table, the 
pressure dependence of T^ in the He-gas experiments is very nearly the same whether 
MgB2 is in single- or polycrystalline form. 
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In Fig. 3 we extend the He-gas results on the same "type B" single crystal to 
much higher pressures using a He- loaded DAC. The initial Tc{P) dependence agrees 
well with the data in Fig. 2, showing a positive curvature at higher pressures. As 
seen in Fig. 3, this curvature is removed if Tc is plotted versus the relative unit-cell 
volume V/Vo The linearity of the Tc versus V/Vo data over such a wide range 

of pressure is remarkable. 

In Fig. 4 are shown our Tc{P) data on an isotopically pure ^^B polycrystalline 



sample to 32 GPa from three separate experiments in a He- loaded DAC [42|; the 



sample used is from the same synthesis batch as in our earlier He-gas studies to 0.7 
GPa • Tc decreases monotonically and reversibly with pressure from ~ 39 K at 



ambient pressure to ~ 15 K at 32 GPa. The initial slope, dTc/dP ~ -1.1 K/GPa, 
is the same as in the He-gas data. As will be discussed below (see Fig. 7), the 
single-crystal and polycrystalline data are in good agreement. 

The only other Tc{P) measurements on polycrystalline MgB2 to very high pres- 



sures using dense He as pressure medium were carried out by Struzhkin et al. 
and are included in Fig. 4 for the ^^B isotope; these authors use a double-modulation 
ac susceptibility technique which determines the superconducting onset rather than 
the superconducting midpoint. As seen in Fig. 4, the agreement with our data 
is remarkably good to 20 GPa, but begins to deviate at higher pressures. Parallel 
studies by the same authors |^ to 44 GPa in dense He on a ^°B isotopic sample yield 
a similar Tc{P) dependence which lies ~ 1 K above their ^^B data below 20 GPa, but 
gradually merges at higher pressures. Unlike our data in Figs. 3 and 4, where Tc is 
seen to be a linear function of V/Vo-, the data of Struzhkin et al. |H3|] show a break in 
slope dTc/ dV near 15-20 GPa. For a full discussion of the latter data see the paper 
by Goncharov and Struzhkin in this special edition. 

Non-He-gas studies to very high pressures include those of Tang et al. to 9 



GPa using a cubic-anvil cell with Fluorinert pressure medium; the Tc{P) dependence 
is in good agreement with that in Fig. 4. However, in DAC studies using methanol- 
ethanol, Tissen et al. found a much larger initial slope {dTc/dP ^ -2 K/GPa) 
accompanied by a relatively large drop in Tc to approximately 6 K at 28 GPa; these 
authors interpret a break in slope dTc/dP near 10 GPa as evidence for a topological 
transition. We note that the present Tc{P) data lie 1-2 K below those of Razavi 
et al. 1^ to 11 GPa obtained using steatite pressure medium. In further DAC 
measurements using steatite Bordet et al. |^ and Monteverde et al. report 
Tc{P) dependences for four different samples which generally lie well above those in 
the present measurement. 
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4 Discussion 



4.1 Intrinsic Dependence of Tc on Pressure 

Before attempting a quantitative analysis of the present data, we would like to first 
discuss possible origins for the differing values of dTc/dP for MgB2 reported in the 



literature (see the Table). Tissen et al. IHSll and, somewhat later, Lorenz et al. H7 



have presented data indicating a strong inverse correlation between the magnitude 
of the initial slope dTc/dP and the value of Tc at ambient pressure, i.e. \dTc/dP\Q 
is larger if Tc(0) is smaller. To reexamine this possible correlation, we have plotted 
in Fig. 5 the initial slope dT^/dP using the high-pressure data in the Table. Taken 
as a whole, the data in Fig. 5 would appear to give some support to the proposed 
strong inverse correlation. Such a correlation is, however, not supported by the 
hydrostatic He-gas data. With the exception of the single data point of Lorenz et al. 
|p9| , all known {dTc/dP)o values for single- or polycrystalline samples obtained using 
He pressure medium lie between -1.07 and -1.2 K/GPa. Although the He-gas data 
do not support a strong correlation between \dTc/dP\Q and Tc(0), samples with lower 
values of Tc(0) (37.16 K versus 39.1 K) do appear to exhibit slightly (~ 5%) larger 
initial slopes dT^/dP (-1.17 K/GPa versus -1.11 K/GPa). In contrast, the {dTc/dP)^ 
values obtained using other less hydrostatic pressure media are often larger (up to 
50-80%) in magnitude. This indicates that differences in the values of {dTc/dP)o 
in the literature may depend more on the pressure medium used than on differences 
between samples, at least as reflected in their Tc(0) values. 

One way to help resolve this issue is to carry out parallel Tc(P) measurements on 
the same sample using different pressure media. We carried out two such studies. In 
the first, we placed a sample provided by V. Tissen in our He-gas system, obtaining 
under purely hydrostatic pressure conditions dT^/dP ~ -1.2 K/GPa [EH]. This value 



is 40% less than that (-2 K/GPa) from Tissen et al. 's DAG study |^5| on the same 



sample with methanol/ethanol pressure medium (see the appropriate vertical arrow 
in Fig. 5). In the second experiment we took the same isotopically pure ^^B sample 
studied previously in both our He-gas and He-loaded DAG (where dTc/dP ~ -1.1 
K/GPa) and replaced the He pressure medium with Fluorinert, obtaining the DAG 
data shown in Fig. 4 with initial slope dTc/dP ~ -1.6 K/GPa, a nearly 50% increase 
in slope (see the appropriate vertical arrow in Fig. 5). Similar effects are observed 



with single crystals; the slope -2.0 K/GPa observed in Fluorinert is much larger 
than that (-1.1 K/GPa) found in the present He-gas study on crystals from the same 
source. The large values of \dTc/dP\Q reported in the literature thus appear to arise 
from the use of fluid pressure media, such as Fluorinert or methanol/ethanol, which 
freeze solid at temperatures well above ~ 40 K. The present results thus do not 
support the existence of the originally proposed [^, ^ strong inverse correlation 
between \dTc/dP\Q and Tc(0). 



In is interesting to note that the DAG measurements of Schlachter et al. |30 
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with solid NaF pressure medium yield a large negative slope with increasing pressure 
(-1.6 K/GPa) which reduces to -1.13 K/GPa as the pressure is reduced to ambient, 
leaving an ambient pressure value Tc(0) permanently suppressed by ~ 12% from the 
initial value. The authors infer that, because of shear stress effects from the solid 
pressure medium, the sample is degraded when pressure is applied, with no further 
degradation upon pressure release. It is also noteworthy that the value dTc/dP ^ 
-1.13 K/GPa is the same to 0.4 GPa as that obtained on the same sample in a He-gas 
experiment by the same group (see Table). 

With these results in mind, it is difficult to understand the relatively small values 
of \dTc/dP\Q reported in the resistivity studies by Monteverde et al. and Bordet 



et al. using solid steatite pressure medium. Since their temperature-dependent 
resistivity data at different pressures was not published, it is difficult to speculate on 
possible origins for this difference. The use by these authors of the resistivity onset 
to define Tc could lead to substantial errors in estimating {dTc/dP)o, particularly if 
significant broadening occurs in the transition under pressure. We also note that in 
this type of pressure cell it is particularly difficult to obtain reliable resistivity data 
and pressure values in the lower 1-2 GPa pressure range. In these papers [0, |46 



no information was given whether Tc(P) was reproducible for identical samples, re- 
versible in pressure, or whether the transition broadened significantly under pressure. 
We note that Struzhkin et al. remeasured Tc{P) in their DAG with no pressure 
medium whatsoever, pressing the diamond anvils directly onto the sample and gasket. 



and reported results resembling the two lower curves of Monteverde et al. It 
is also noteworthy that Razavi et al. also used steatite pressure medium, and 
obtained the initial slope -1.03 K/GPa; in their published resistivity data there is 
little change in the transition width with pressure so that the estimate of the shift in 
Tc is relatively straightforward. 

In a further experiment, we left out the pressure medium entirely and pressed two 
WG Bridgman anvils with 6 mm dia. flats for 5 minutes directly onto ~ 5 mg MgB2 
powder with 10 tons force; this typically results in a pressure distribution across the 
center of the disc resembling a bell-shaped curve with P ^ 8 GPa at the center 
and P = at the edge. The sample was thoroughly compacted by this procedure, 
resulting in a disc-shaped sample approximately 70 fim thick. After removing the 
compacted disc from between the anvils, the disc was separated into three separate 
regions (center, middle, and outside), as seen in Fig. 6; the sample from each region 
was then gently broken up and placed into our ac susceptibility coil system. As 
seen in Fig. 6, the sharp superconducting transition of the virgin powder sample 
is broadened by a significant amount (~ 10 K) by the compaction procedure, the 
broadening being somewhat less from the center to the middle to the outside. This 
broadening is presumably the result of strong plastic deformation and the resulting 
internal strains in the material. If strong internal strains can lead to broadening as 
large as 10 K at ambient pressure, it is not unreasonable to assume that such strains 
are capable of causing enhanced downward shifts in under nonhydrostatic pressure 
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conditions. 

The important role that defects and strains play in MgB2's superconducting state 
is emphasized by the fact that the value of Tc in a large number of thin-film and bulk 
MgB2 samples appears to follow a Testardi correlation [0, Tc is lower for samples 



in which the conduction electrons are strongly scattered. Lorenz et al. also report 
that the value of is degraded with increasing lattice strain, a result confirmed by 
Serquis et al. . in MgB2 has also been found to be lowered following mechanical 



milling [50| and after irradiation by fast neutrons |51|. In contrast to the work of 



Lorenz et al. and Serquis et al. , Hinks et al. have recently reported that 
accidental impurity doping can have an effect on Tc much larger than that of lattice 
strain. Additionally, they show that grain-interaction stresses can significantly alter 
the lattice parameters of MgB2 depending on the impurity phases present with the 
MgB2. For example, the largest strains are seen for samples where MgB4 is present 
as an impurity phase, imparting to the bulk sample the properties of an MgB2/MgB4 
composite. One might hypothesize that grain-interaction stresses could significantly 
alter the response of MgB2 grains to pressure in such samples, and, in this way, modify 
the observed response of to pressure. 

In contrast to the poor agreement between the large body of nonhydrostatic Tc{P) 
data on MgB2, the agreement between the purely hydrostatic He-gas data to 0.7 GPa 
and the nearly hydrostatic dense He data to 20 GPa (see Fig. 4) is remarkably 
good for both single-crystalline and polycrystalline samples. It is thus reasonable 
to assert that the intrinsic initial pressure dependence for MgB2 samples with the 
highest values of Tc (39.1 K for ^^B) is given by dTjdP ~ -1.11(2) K/GPa, the T^{P) 
dependence to 32 GPa being given by the present data shown in Figs. 3 and 4. We 
would now like to compare these experimental results with theoretical models. 

4.2 Comparison with Theory 

For most known superconductors, including MgB2, Tc is found to decrease with pres- 
sure. In fact, dTjdP is negative for all simple-metal superconductors (e.g. Pb, Al, 



Sn, and In [^) due to pressure-induced lattice stiffening (higher phonon frequencies) 



which weakens the electron-phonon coupling RbaCeo and other alkali-doped 



fullerenes are exceptions to this "rule" : here the rapid decrease in Tc with pressure 
originates not from lattice stiffening but rather from a sharp decrease in the electronic 
density of states N{Ef) |Q. What is the origin of the negative pressure dependence 
of Tc for MgB2 - is it lattice stiffening, a decrease in the electronic density of states, 
or something else? 

Theoretical models are calculated in terms of the dependence of the relevant prop- 
erties on the lattice parameters or unit cell volume V. From the intrinsic initial slope 
dTc/dP ~ -1.11 K/GPa for MgB2 we can calculate the logarithmic volume derivative 



dlnTr B f dT, 



d\nV Tc(0) V^^^ 



+4.18, (1) 
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where we use Tc(0) = 39.1 K from above and the bulk modulus B = 147.2 GPa from 
He-gas neutron diffraction studies on the same sample 

We now convert the Tc{P) data for polycrystalline MgB2 in Fig. 4 to Tc{V) data 
using the Murnaghan equation of state In Fig. 7 we compare the resulting 



Tc{V) dependences to that for the type "B" single crystal from Fig. 3; the agreement 
is excellent considering that Tc(0) for the single crystal is nearly 1 K lower than for 
the polycrystalline sample. 

We now attempt a fit to the data in Fig. 7 using the well known McMillan formula 



T^Mexp/ -1-04(1 + A) 1 
1.2''''PlA-^*(l + 0.62A)/' 

which is valid for strong coupling (A < 1.5) and connects the value of with the 
electro n-phonon coupling parameter A, an average phonon frequency (u) , and the 
Coulomb repulsion fi*. We should not expect too much from this fit since the 
McMillan formula is a solution of the isotropic Eliashberg equations and thus ignores 
the strong anisotropics in the vibrational, electronic, and superconducting properties 
of MgB2. In the following we draw heavily on the detailed analysis of Chen et al. 



| 55| who consider the effect of pressure on the three relevant parameters in Eq. (2) 
{oj) , A, and /x*, . 

Taking the logarithmic volume derivative of both sides of Eq. (2), we obtain the 
simple relation 

-7-Ai in7T7 +^2 7tt:T7 + 27 , (3) 



d\nV ' \ d\nV J {dlnV 

where 7 = —d\n{(jj) /d\n.V is the Griineisen parameter, rj = N{Ef) (P) is the 
McMillan-Hopfield parameter given by the product of the electronic density of states 
and the average squared electronic matrix element, and the dimensionless prefac- 
tors are given by Ai = 1.04/i*(l + A) [1 + 0.62A] / [A - /i*(l + 0.62A)]^ and A2 = 
1.04A [1 + 0.38/i*] / [A - /i*(l + 0.62A)]^ . Both Ai and A2 are calculated using the 
values of A and /i* at ambient pressure. 

Eq. (3) looks rather formidable, but has a very simple interpretation. Fortunately, 
the first two terms on the right are usually small relative to the third term, as we 
will see below, so that dlnTc/dlnV ~ A2{dlnr]/d\nV + 27}. Since A2 is always 
positive, the sign of the logarithmic derivative d\nTc/d\nV is determined by the 
relative magnitude of the two terms dlnrj/dlnV and 27. The first "electronic" term 
is negative {dlnrj/dlnV ~ —1 for simple metals (s,p electrons) |2^, but may equal 
-3 to -4 for transition metals (d electrons) ||5^), whereas the second "lattice" term is 
positive (typically 27 ^ 3 — 5). Since in simple-metal superconductors, like Al, In, Sn, 
and Pb, the lattice term dominates over the electronic term, the sign of dlnTc/dlnV 
is the same as that of {dlnrj/dlnV + 27}, namely positive; this accounts for the 
universal decrease of Tc with pressure due to lattice stiffening in simple metals. In 
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selected transition metals, the electronic term may become larger than the lattice 
term, in which case would be expected to increase with pressure, as observed in 
experiment p6 |. 



We now apply the McMillan equation to the above results on MgB2 using at 
ambient pressure the logarithmically averaged phonon energy from inelastic neutron 



studies [0 {lu) = 670 K, T^{0) ~ 39.1 K, and /x* = 0.1 |8|, yielding A ~ 0.898, Ai = 
0.558 and A2 = 1.76. From the expression derived by Chen et al. for s,p metals 
d\nT]/d\nV = -[d\nN{Ef)/d\nV] - 2/3 and the value d\nN{Ef)/d\nV = +0.46 
from Loa and Syassen |Q , the dependence of the Hopfield parameter is estimated to 
be d\nT]/d\nV = —1.13, a value reasonably close to the generic value (-1) used above 
in the analysis for the simple-metal superconductors and to the value (-0.81) obtained 
from first-principles electronic structure calculations on MgB2 by Medvedera et al. 
||60|| . Chen et al. find that /i* increases only weakly with pressure at a rate 
(f) = din fi* /dlnV = -O.I7— 0.035; note that this derivative is very small (~ -0.3) so 
that the second term on the right side of Eq. (3) is relatively unimportant. We now 
have estimates of all quantities on the right side of Eq. (3) except 7 which we use 
as a fit parameter. Setting the left side of Eq. (3) equal to the experimental value 
-1-4.18 from Eq. (1), we find 7 = 2.39, in reasonable agreement with the value 7 ~ 2.9 



from Raman spectroscopy studies |£T[ or 7 ?a 2.3 from ab initio electronic structure 
calculations on MgB2 |6^. Note that for the present "type B" crystal Tc(0) ~ 38.24 
K and dTjdP ^ -1.10 K/GPa which gives dlnTjdlnV = +4.23, A ~ 0.887, and 
7 = 2.36. 

Eq. (3) is only valid for small changes in the parameters, i.e. for experiments to 
a few GPa pressure where the change in unit cell volume is only a few percent. To 
attempt to fit the very high pressure Tc{V) data in Figs. 3 and 7, we need to use 
the full McMillan equation and insert explicitly the change in the parameters with 



relative volume. As suggested by Chen et al. [^, we set 

{cu) = {co), iV/Vo)-\ A = XoiV/Vor, and /x* = fi*oiV/Vo)^ (4) 

where cp = {dlnr]/dlnV + 27}. Using the values of the parameters for the sample 
with Tc(0) = 39.1 K used above to fit the initial pressure dependence (7 ~ 2.4, 
(f = -1.13 + 2 X 2.4 = 3.67, and (f) = -0.1 x 2.4 - 0.035 = -0.275 and inserting 
the volume dependences from Eq. (4) into the McMillan equation, the lower solid 
fit curve in Fig. 7 is obtained. Note that this curve clearly lies below the data at 
higher pressures. A reasonably good fit to the data over the pressure range to 32 
GPa {V/Vo = 0.855) is found for 7 = 2.2 (see Fig. 7). As with the simple s,p 
metal superconductors, Tc in MgB2 appears to decrease under pressure due to lattice 
stiffening. Note that within this approximation approaches K asymptotically 
at very high pressures. The fit curve for 7 = 2.2 is predicted to fall below 1 K for 
V/Vo = 0.73 which corresponds to an applied pressure of 93 GPa. The fact that the 
experimental data can be well fit by the McMillan formula with reasonable values of 
the parameters lends support to the view that MgB2 is a BCS superconductor with 
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moderately strong electron-phonon coupling. However, the relatively small value of 
the fit parameter 7 = 2.2 compared to experiment (2.9 ||6l|) is cause for concern. 

As discussed in the Introduction, the binary compound MgB2 is a quasi 2D system 
with highly anisotropic electronic and lattice properties, including multiple supercon- 
ducting gaps. The above analysis of experimental data using the McMillan formula, 
which represents a solution to the isotropic Eliashberg equations, is a good first step 
but does not permit unequivocal conclusions. What is needed is an extension of 
the solution of the fully anisotropic Eliashberg equations [|13| to reduced lattice pa- 
rameters. The initial dependences dTc/dP ~ —1.11(2) K/GPa and dlnTc/dlnV ~ 
-1-4.18, and the Tc{P) and TdV) data to 32 GPa in Figs. 3, 4, and 7 stand ready to 
provide a stringent test of such a calculation. 
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5 Figure Captions 



Fig. 1. Real part of the ac susceptibility x' (1023 Hz, 1 Oe (rms)) versus tempera- 
ture for a "type B" MgB2 single crystal at three different pressures. The position of 
the superconducting midpoint is indicated for ambient pressure data. 
Fig. 2. Dependence of the superconducting transition temperature of one "type 
A" (•) and two "type B" (•, x) MgB2 single crystals on hydrostatic He-gas pres- 
sure. Numbers give order of measurement. Pressure was normally changed at 
room temperature; for data with primed numbers the pressure was changed at low 
temperatures (~ 50 K). Solid and dashed lines are guides to the eye. 
Fig. 3. Dependence of the superconducting transition temperature of "type B" 
MgB2 single crystal on nearly hydrostatic pressure in a He-loaded DAC. All data 
taken in order of increasing pressure, where pressure was only changed at ambient 
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temperature. is determined from the midpoint of the superconducting transition 
in the temperature-dependent ac susceptibihty x' (1003 Hz, 3 Oe (rms)). Dashed hne 
gives slope dTc/dP ~ —1.10 K/GPa of He-gas measurement on same crystal (Fig. 2). 
The "error bars" give the temperatures of the onset and end of the superconducting 
transition; note that the transition broadens for P > 14 GPa. Tc versus relative 
volume V/Vo is also shown; straight sohd line is guide to the eye. 
Fig. 4. Dependence of the superconducting transition temperature of isotopically 
pure ^^B polycrystalline MgB2 on nearly hydrostatic pressure in He-loaded DAC 
measurements: three different experiments (•,■,□, A), closed symbols (increasing 
pressure), open symbols (decreasing pressure); measurements from Struzhkin et al. 
p3| (x). Measurements in our DAC on same sample with Fluorinert pressure medium 
(♦) are also shown; solid and dashed lines gives slopes (-1.11 K/GPa) and (-1.6 
K/GPa) from He-gas |27| and Fluorinert data, respectively, on same sample. 
Fig. 5. Initial pressure derivative {dTc/dP)o versus value of Tc at ambient pressure; 
data are taken directly from the Table. Measurements with He (•, o) and non- 
He (-I-, x) pressure media; both polycrystalline (•, +) and single- crystalline (o, x) 
samples are represented. Vertical arrows show change in measured value of {dTj dP)o 
for a given sample upon changing to He pressure medium (see text). 
Fig. 6. Relative change in the real part of the ac susceptibility x' of MgB2 versus 
temperature for both loose powder and flat, compacted samples. For the compacted 
disc, samples were taken from the center, middle, and outside regions. 
Fig. 7. Tc values from present single-crystalline (open stars) and polycrystalline 
(•, ■, □, a) data from Figs. 3 and 4 plotted versus relative volume. Solid lines are 
calculated curves using McMillan's equation (Eq. 2) for three different values of the 
Griineisen parameter 7 (see text). 
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Table. Summary of available high-pressure Tc{P) data on MgB2 single-crystals (first 
4 rows) and polycrystals (remaining rows). T^. values are at ambient pressure from 
the superconducting midpoint in the ac susceptibility Xac and electrical resistivity p 
measurements; Struzhkin et al. |Q use a double-modulation technique xSf'^ which 
is beheved to give the superconducting onset. {dTc/dP)o is the initial pressure 
derivative. P™'^^(GPa) is the maximum pressure reached in the experiment . Unless 
otherwise specified, samples with the natural boron isotopic abundance lO-^^B are 
studied. Arrows indicate increasing (t) or decreasing (|) pressure. 
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